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Abstract Image analysis was used in this study to develop a 
faster, more objective method for quantitatively measuring 
tracheid shapes on cross section. The morphological fea¬ 
tures were measured automatically using the circularity in¬ 
dex and aspect ratio based on the skeletonized images. The 
rectangular and hexagonal similarities were introduced to 
evaluate the tracheid shapes. The results for 18 species of 
Japanese softwoods showed that the difference in tracheid 
shapes among the species could be successfully estimated 
using the quantitative method. In addition, the boundary of 
latewood might be evaluated based on the change in circu¬ 
larity indexes. 

Key words Tracheid shape ■ Circularity index ■ Aspect 
ratio ■ Shape similarity • Image analysis 


Introduction 

Image-processing techniques have been widely used in 
quantitative analyses of cells or particle shapes in various 
fields, for example, automatic classification of crops accord¬ 
ing to morphology, 1 determining the optimum operating 
condition of the impactor particles of cement powder, 2 and 
predicting the concentration of airborne dust in woodwork¬ 
ing chambers. 3 For studying wood anatomy, measurements 
of cell shapes have become possible. Gonzalez and Wintz 4 
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and Levine 5 noted that several stand edge detectors are 
available for delineation of cell boundaries using computer 
image-processing techniques. As a new approach, image 
processing based on the distance transform 6 was used for 
quantitative assessment of wood cell boundaries. Recently, 
element angles (angles between two radial walls) of trache- 
ids and variation of shape in a single annual ring were 
investigated using fast Fourier transform (FFT) image 
analysis. 7,8 

In our previous study 9 the morphological features of 
cross-sectional vessel shapes for 47 Japanese species were 
analyzed quantitatively using image processing. The circu¬ 
larity index has proven to be a useful means for assessing 
cell shape difference between solitary pores and pore mul¬ 
tiples. The variability of cross-sectional tracheid shapes, 
especially earlywood, is not large; but recently some studies 
have shown that there is an important effect of cross- 
sectional tracheid shapes on the anisotropy of wood proper¬ 
ties. 1011 Therefore, studying the morphological features of 
tracheids is necessary to explain the anisotropy or to predict 
various properties of wood using theoretical models. In this 
study we characterized the morphological features of trac¬ 
heid shapes on cross sections and their changes within an 
annual ring using a circularity index and an aspect ratio by 
an image-processing technique. 


Materials and methods 

Microscopy 

Samples from 18 Japanese softwoods were used in this 
study, as shown in Table 1. Measurements of morphological 
features were made by different annual ring widths. Cross 
sections, 10-20 pm thick, were cut with a sliding microtome, 
stained with 1% aqueous safranine, and mounted on glass 
slides. Photographs were taken continuously along the ra¬ 
dial direction from earlywood to latewood in every annual 
ring with a 20X objective and an FK2.5X ocular. The mag¬ 
nification of the photographs was 364. 
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Table 1 . Morphological features of tracheids in earlywoods for 18 Japanese softwoods 


Species 

Scientific names 

RW (mm) 

LP (%) 

C/ AV 

2 A v 

S T 

S H 

No. 

Monti 3 

Abies firma Sieb. et Zucc. 

1.27 

14.60 

0.782 (0.094) 

1.532 (0.746) 

94.09 (4.70) 

92.10 (6.59) 

3930 

Todomatsu 

Abies sachalinensis Mast. 

1.17 

10.05 

0.728 (0.157) 

1.708 (0.699) 

92.18 (9.24) 

89.66 (11.45) 

2088 

Hinoki 

Chamaecyparis obtusa Endl. 

4.03 

2.46 

0.822 (0.119) 

1.733 (0.502) 

88.48 (5.10) 

94.57 (7.59) 

8703 

Sawara 

Chamaecyparis pisifera Endl. 

1.19 

6.05 

0.819 (0.091) 

1.753 (0.848) 

90.47 (3.51) 

96.09 (4.44) 

3761 

Sugi 

Cryptomeria japonica D. Don 

1.97 

7.48 

0.829 (0.087) 

1.638 (0.744) 

90.82 (3.02) 

95.23 (3.08) 

10994 

Karamatsu 

Larix kaempferi Carr. 

1.79 

15.65 

0.821 (0.081) 

1.613 (0.469) 

91.15 (3.66) 

96.41 (3.97) 

3486 

Ezomatsu 

Picea jezoensis Carr. 

1.25 

16.36 

0.826 (0.130) 

1.582 (0.724) 

87.93 (6.10) 

93.92 (8.91) 

2868 

Akamatsu 3 

Pinus densiflora Sieb. et Zucc. 

1.75 

25.03 

0.841 (0.075) 

1.566 (0.501) 

89.96 (3.68) 

94.87 (3.10) 

3880 

Himekomatsu 

Pinus parviflora Sieb. et Zucc. 

1.14 

19.16 

0.797 (0.093) 

1.650 (1.014) 

93.64 (2.50) 

94.25 (3.27) 

2747 

Kuromatsu 

Pinus thunbergii Pari. 

1.01 

18.57 

0.822 (0.084) 

1.598 (0.435) 

91.99 (3.26) 

93.80 (3.07) 

8700 

Inumaki 

Podocarpus macrophylla 

D. Don 

1.82 

11.82 

0.777 (0.073) 

1.444 (0.425) 

95.84 (3.57) 

90.54 (4.91) 

8458 

Togasawara 

Pseudotsuga japonica Beissn. 

3.83 

24.33 

0.803 (0.050) 

1.623 (0.081) 

93.02 (2.85) 

94.74 (1.24) 

3771 

Koyamaki 

Sciadopitys verticillata Sieb. 
et Zucc. 

0.75 

14.81 

0.767 (0.072) 

1.518 (0.561) 

96.68 (3.04) 

90.81 (4.41) 

3020 

Ichii 

Taxus cuspidata Sieb. et Zucc. 

0.73 

17.16 

0.817 (0.076) 

1.691 (0.562) 

91.11 (3.26) 

96.52 (3.01) 

3297 

Nezuko 

Thuja standishii Carr. 

1.37 

23.05 

0.819 (0.104) 

1.705 (0.783) 

89.55 (4.71) 

95.30 (5.71) 

3793 

Hiba 

Thujopsis dolabrata Sieb. 
et Zucc. 

1.03 

10.21 

0.821 (0.078) 

1.549 (0.509) 

91.88 (3.31) 

95.95 (3.41) 

2913 

Kaya 

Torreya nucifera Sieb. et Zucc. 

0.75 

20.96 

0.763 (0.072) 

1.513 (0.567) 

94.62 (3.14) 

90.31 (4.21) 

2385 

Tsuga 

Tsuga sieboldii Carr. 

1.31 

25.77 

0.776 (0.066) 

1.509 (0.454) 

96.26 (3.10) 

91.30 (4.13) 

3160 


The data in parentheses are standard deviations. 

RW, average value of five annual ring widths; LP, average latewood percentages for five annual rings; C/ AV , average circularity indexes; 2 AV , 
average aspect ratios; S T , rectangular similarities; S H , hexagonal similarities; No, number of measured tracheids. 

“The number of annual rings is 11. 


Image analysis 

The photographs were digitized by an image scanner (flat¬ 
bed type, IX-4015; Canon) with an image resolution of 64 
dpi. A Macintosh 8100/80AV personal computer system 
was used for storing and analyzing the digital images. When 
input, the horizontal coordinate in a cathode-ray tube 
(CRT) display corresponds to the tangential direction of the 
original image. The monochromatic images with a good 
contrast were 400 X 600 pixels and 256 gray levels. The 
magnification of 364 for the photographs corresponded to 
1.09 pm/pixel of the original images. The morphological fea¬ 
tures of individual tracheids in the digitized images were 
measured using the software NIH (National Institutes of 
Health, USA) Image/ppc 1.56b77. 

Because copylilms were used during photographing, 
photographs with a clear boundary between the cell wall 
and cell lumen were obtained. Monochromatic images with 
good contrast were easily binarized. The binary cell wall 
image produced was then spatially filtered to remove high- 
frequency noises. An opening process was used to remove 
small foreground objects and irregularities in the larger 
object boundaries. Then a closing process of one pixel was 
used to remove small holes from the remaining foreground 
objects (cell wall). 

Although the cell lumina can be clearly distinguished, 
the boundaries between cells are less obvious. To discern 
these boundaries on the assumption that the cell boundaries 
are midway between adjacent cells, we used a shrinkage 
technique that is similar to the distance transform. A skel¬ 
eton of the binary cell wall image with a single pixel ap¬ 
proximates the position of the middle lamella between 


adjacent cell walls. It is possible to measure various param¬ 
eters of the cells from these images, such as the cross- 
sectional area and perimeter. An exception occurs when 
making measurements of cells near the ray cells. When 
shrinkage processing was used, the shrinkage of the ray cells 
caused deviation of cell shapes near ray cells. This was due 
to a precondition of shrinkage processing that lay on the 
assumption that the double wall thickness was the same. 
Thus, measurements of the cell lumen were carried out near 
the ray cells without the shrinkage processing. 


Feature extraction 


Morphological features of particles can be described by 
many parameters, such as area, perimeter, Heywood diam¬ 
eter, maximum length, maximum breadth, aspect ratio, and 
circularity index as a shape factor. In this study, the circular¬ 
ity indexes and aspect ratios as the important morphological 
factors were measured to evaluate quantitatively the shapes 
of tracheids. The circularity index is defined by the follow¬ 
ing formula. 12 


CI(X) 


4%A 

~cF 


( 1 ) 


where Cl is the circularity index, C is the perimeter, A is the 
area, and A is the aspect ratio expressed by the ratio of 
Feret’s radial diameter (Z)jj) and Feret’s tangential diam¬ 
eter (£>t). The Cl is 1.0 for any circle and 0.785 for any 
perfect square, and their aspect ratios are 1.0. 

Ohgoshi et al. 13 used the relation between the circularity 
index and the aspect ratio to evaluate the shapes of 
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Fig. 1. Relation between the circularity indexes of rectangles or hexa¬ 
gons and the aspect ratios. Cl, circularity index; A, aspect ratio; C/(A), 
circularity index of a tracheid; C/ H (A), circularity index of a hexagon; 
C/ T (A), circularity index of a rectangle 



Fig. 2. Relation between the circularity index (Cl) and aspect ratio (A) 
for various models. S T , rectangular similarity; S H , hexagonal similarity 


bordered-pit membrane pores of softwoods by the image 
analysis system. Fujiwara 1415 also used the circularity index 
and aspect ratio (D^JDj) to analyze the changes in fiber 
shapes during fiber differentiation in hardwoods. Usually 
when an effect of tracheid shapes on cross sections on the 
wood properties is discussed, rectangular or hexagonal 
models are used. It is necessary to evaluate whether the 
shapes of tracheids in any species are more similar to a 
rectangle or a regular hexagon. The circularity index of a 
tracheid is compared with that of a rectangle or a hexagon 
with the aspect ratio (A) of a tracheid. A similarity to a 



Fig. 3. Comparison of circularity indexes (Cl) between the lumen and 
middle lamella boundaries. Solid circles, lumen boundary; open circles, 
middle lamella boundary; solid triangles, differences in circularity in¬ 
dexes between lumen and middle lamella boundaries. A Ichii (Taxus 
cuspidata Sieb. et Zucc.). Ring width is 0.45 mm and latewood percent¬ 
age 15.35%. B Sugi (Crytomeria japonica D. Don). Ring width is 
1.87 mm and latewood percentage 18.96%. Dotted lines, boundary of 
latewood according to Mork’s definition 


rectangle or a regular hexagon was introduced into evalua¬ 
tions of tracheid shapes. The rectangular similarity (S T ) and 
hexagonal similarity (S H ) were calculated by the following 
formula: 


[l - |c/(A) - C/ t (A)|] X 100 

(2) 

[l - |C/(A) - C/ h (A)|] x 100 

(3) 


From this formula it is understood that if the circularity 
index [C7(A)] of the tracheid is equal to that of a rectangular 
model [C7 X (A)] or a regular hexagonal model [C/ H (A)], S r or 
S H is 100. This shows that the shape of a tracheid is the same 
as a rectangle or a regular hexagon. 

The circularity indexes of a rectangle [C/ T (A)] and a 
regular hexagon [C/ H (A)] were calculated based on the as¬ 
pect ratio (A) of a tracheid using the following formula: 

Rectangle: C7 X (A) = -—-— (4) 

A + 2 + — 

X 
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0.0 0.2 0.4 0.6 0.8 1.0 


Relative position in one annual ring 
B (Tsuga) 

Fig. 4. Change of circularity index (Cl), aspect ratio (A), and shape 
similarity (S) in tsuga (Tsuga sieboldii Carr, left) and momi (Abies 
firma Sieb. et Zucc, right). Solid circles, circularity index; open circles. 


A1 + VI H~ 3A 2 

9 = 2 arcsin—-——— (5) 

Hexagon: < 2[1 + A J 

= -^2sin-^ + sinSj (6) 

where 9 is an element angle 16 (the angle between two radial 
walls) obtained from the aspect ratio (Eq. 5). Figure 1 
shows the relations between the circularity indexes of rect¬ 
angles or hexagons and the aspect ratios. 

To help understand the shape similarity concept, the 
various possible models of earlywood tracheids are shown 
in Fig. 2. Models A, B, and D were assumed to be tracheids 
near rays. Models C and E were assumed to be tracheids in 
earlywood. Their rectangular and hexagonal similarities 
were estimated as shown in Fig. 2. The results show that the 
shapes of the models can be evaluated quantitatively on the 
basis of the shape similarity. 

In addition, when the areas, perimeters, and aspect ratios 
of tracheids were measured, the coordinates of tracheids in 



0.0 0.2 0.4 0.6 0.8 1.0 



0.0 0.2 0.4 0.6 0.8 1.0 

Relative position in one annual ring 
(Momi) 


aspect ratio; solid squares, rectangular similarity; open triangles, hex¬ 
agonal similarity; dotted lines, boundary of latewood according to 
Mork’s definition 


a computer display were also measured. Therefore, one 
annual ring was divided into 10 parts on the basis of the 
coordinates from earlywood to latewood to obtain the dis¬ 
tribution of circularity indexes and the shape similarities. 

Accuracy of a circularity index 

In a previous paper 9 we made circles with various diameters 
of about 5-37 pixels using a paint software to discuss the 
effects of digitizing processing on the accuracy of a circular¬ 
ity index. According to the results obtained, the circularity 
index of the circle increased rapidly as the circle diameter 
increased. If the diameter of a circle was larger than 35 
pixels, the relative error of the circle perimeters was less 
than 3%, and that of the circularity indexes of the circles 
was less than 5%. To reduce the errors that resulted 
fr om the digitizing process, the Heywood diameters (d = 
\ / 4Aln) of tracheids were estimated. The circularity in¬ 
dexes of tracheids were corrected on the basis of the 
Heywood diameters using the regression formula of a mea¬ 
sured circularity index and its diameter (Cl = 1.039 - 0.315/ 
lnrf). 
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Relative position in one annual ring 
B (Nezuko) 


Relative position in one annual ring 
(Akamatsu) 


Fig. 5. Change of circularity index (Cl), aspect ratio (2), and shape similarity (5) in nezuko (Thuja standishii Carr, left) and akamatsu (Pinus 
densiflora Sieb. et Zucc, right). Symbols are the same as in Fig, 4 


To determine the effects of a threshold gray level during 
binarizing on the circularity indexes, the average values of 
about 100 earlywood tracheids were measured in the differ¬ 
ent threshold gray levels for three species: sugi ( Cryptom- 
eria japonica D. Don), todomatsu ( Abies sachalinensis 
Mast.), and akamatsu ( Pinus densiflora Sieb. et Zucc.). 
There are few variations of the circularity indexes from 
threshold values of 100-200. This is due to monochromatic 
photographs having a good contract when a copyfilm was 
used. For all monochromatic images, the threshold gray 
level was about 150 with binary processing. 


Results and discussion 

Shapes of middle lamella boundary and lumen 

To compare the difference of shapes between the middle 
lamella boundary (MLB) and the lumen, their circularity 
indexes were measured. Figure 3 shows the shape variations 
of the MLB and lumen for ichii ( Taxus cuspidata Sieb. et 
Zucc.) and sugi ( Cryptomeria japonica D. Don). The results 
show the larger differences in the circularity indexes of the 


MLB and lumen for ichii and the small differences for sugi. 
Because the cell wall thickness in the earlywoods is small, 
the effects of cell corners are relatively small. Generally, as 
the cell walls thicken, the effects of cell corners be¬ 
come larger. On the other hand, tracheid diameters may 
affect tracheid shape. Peachey and Osborne 6 indicated that 
the difference in shapes between the MLB and lumen 
appeared as an increase in cell wall thickness. The cross- 
sectional shapes of tracheids are the shapes of MLB in this 
study. 

Averages of circularity index, aspect ratio, 
and shape similarities 

It was reported that the shapes of tracheids have a larger 
effect in earlywood than in latewood in terms of the pro¬ 
perties of the wood, such as shrinkage and Young’s modu¬ 
lus; the effects of the cell wall structure itself become 
greater as the cell walls thicken. 10,16 Table 1 shows the 
morphological features of tracheids, such as the average 
circularity indexes and aspect ratios and the rectangular 
and hexagonal similarities measured in earlywood of five 
annual ring widths. For akamatsu ( Pinus densiflora Sieb. 
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Fig. 6. Microphotographs of transverse sections of momi (Abiesfirma Sieb. et Zucc) (A) and akamatsu (Pinus densiflora Sieb. et Zucc.) (B). Bars 
200 pm 


et Zucc.) and todomatsu (Abies sachalinesis Mast.), the av¬ 
erages of the circularity indexes are 0.841 and 0.728; 
these are maximum and minimum values, respectively, 
compared with those of other species. For all species, the 
averages of the aspect ratios vary from 1.50 to 1.75. Based 
on the standard deviations, the variations in aspect ratios 
are large. 

The rectangular and hexagonal similarities for each trac- 
heid were calculated according to Eqs. (2) and (3). The 
averages and standard deviations are shown in Table 1 for 
18 Japanese softwood species. Because there is a difference 
between the rectangular and hexagonal similarities, it is 
suggested that similarity is a useful tool for evaluating the 
shape of cross-sectional tracheids. Comparing their shape 
similarities, the rectangular similarities are larger than the 
hexagonal ones for momi (Abies firma Sieb. et Zucc.), 
inumaki (Podocarpus macrophyllus D. Don), koyamaki 
(Sciadopitys verticillata Sieb. et Zucc.), kaya (Torreya 
nucifera Sieb. et Zucc.), and tsuga ( Tsuga sieboldii Carr.). 
The other species have larger hexagonal than rectangular 
similarities. The differences in tracheid shapes among spe¬ 
cies could be evaluated from the results of the two types of 
shape similarity. It is important to explain the production of 
wood anisotropy and to predict various properties of wood 
using theoretical models. 


Change in circularity index and shape similarity within an 
annual ring 

Generally, the variation in tracheid shapes can be examined 
by visual observation, and it is usually believed that tracheid 
shapes are similar to hexagons in early wood and to rect¬ 
angles in latewood. To evaluate this variation quantita¬ 
tively, we measured the changes in circularity index and 
shape similarity within an annual ring. Typical results are 
shown in Figs. 4 and 5. Figure 4B shows that the rectangular 
similarities are larger than the hexagonal similarities in 
tsuga (Tsuga sieboldii Carr.) and momi (Abiesfirma Sieb. et 
Zucc.). It is concluded that the tracheid shapes of the two 
species are closer to rectangles. In contrast, the tracheid 
shapes are closer to hexagons for nezuko (Thuja standishii 
Carr.) and akamatsu (Pinus densiflora Sieb. et Zucc.) (Fig. 
5B). The shapes of many tracheids may be similar to models 
A, B, and D in Fig. 2 because the circularity index of momi 
is close to that of a rectangle (Fig. 4B). Figures 4A and 5A 
demonstrated that there are larger circularity indexes in the 
transition zone than those in earlywood and latewood. 
Jagels and Dyer 17 studied the tracheid shape change of red 
spruce within an annual ring. Their results showed that the 
boundary between earlywood and latewood could be evalu¬ 
ated based on the peak of the circularity index. In Figs. 4 
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Fig. 7. Relation between shape similarities (S) and annual ring width Open squares , rectangular similarity; open circles : hexagonal similarity, 
as well as the latewood percentage for momi (Abies flrma Sieb. et Vertical lines represent standard deviations. R T and R H , correlation 
Zucc.) (left) and akamatsu (Pinus densiflora Sieb. et Zucc.) (right), coefficients 


and 5, the dotted lines show the boundary of latewood 
obtained according to Mork’s definition. 18 ' 19 There is a good 
agreement between the boundary of latewood based on 
Mork’s definition and the peak of the circularity index 
rather than the shape similarities. 

Effects of annual ring width and latewood percentage 

To determine the effects of annual ring widths and late¬ 
wood percentages on the similarities of tracheids, we mea¬ 
sured the shape similarities for momi (Abies firma Sieb. et 
Zucc.) and akamatsu (Pinus densiflora Sieb. et Zucc.) in 11 
annual ring widths. According to the visual observation 
(Fig. 6), there are clear differences in tracheid shape be¬ 
tween the two species. 

Figure 7 shows the relation between the shape similarity 
and annual ring width and the latewood percentage for 
momi and akamatsu. The shape similarity, as the average of 
all tracheids in every ring and its standard deviation, were 
obtained. The annual ring width and latewood percentage 
were measured, when measuring the shape similarity of 


tracheids in every ring. There was little correlation between 
the rectangular similarity and annual ring width because 
there are smaller correlation coefficients at the 0.95 con¬ 
fidence level. There is a higher correlation between rectan¬ 
gular similarity and latewood percentage for momi. Flere 
the contribution of the latewood tracheids with the larger 
rectangular similarities become larger as the latewood per¬ 
centages increase. With this increase, it is evident that the 
number of latewood tracheids with higher rectangular simi¬ 
larity increases. The results for akamatsu show that the 
rectangular similarity and hexagonal similarity have little 
change as the latewood percentage increases. This is be¬ 
cause the shapes of latewood tracheids were similar to hexa¬ 
gons, although the latewood percentage increased. In 
addition, based on the results for akamatsu in Fig. 5, it was 
proved that the tracheid shapes in latewood are similar to 
hexagons because the hexagonal similarities in latewood 
have large values. For other species, the effects of the an¬ 
nual ring width and the latewood on the shape similarity 
cannot be discussed because of few number of their annual 
rings. 
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Conclusion 

Using image analysis techniques, the morphological proper¬ 
ties of tracheids on cross section for 18 species of Japanese 
softwoods were measured quantitatively based on the circu¬ 
larity index and aspect ratio. The rectangular and hexagonal 
similarities were introduced to evaluate the tracheid shapes. 
The difference in tracheid shapes among species could be 
successfully estimated using this quantitative method. In 
addition, the boundary of latewood might be evaluated us¬ 
ing the change in circularity index. However, using the 
shape similarity concept, it would be difficult to reconstruct 
the models of tracheid shapes because the cell wall dimen¬ 
sions (e.g., cell wall length and the element angle) as impor¬ 
tant parameters of tracheid shape could not be determined. 
It is necessary to measure automatically the dimensions and 
element angles of cell walls using the FFT method. 
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